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The potentials between static sources in various representations in SU(3) are calculated based 
on the fat-center-vortices model of Faber, Greensite and Olejnik. At intermediate distances, most 
distributions of the flux within vortices lead to potentials that are qualitatively in agreement with 
"Casimir scaling," which says that the string tension is proportional to the quadratic operator of 
the representation. However, at the quantitative level, violations of Casimir scaling are generally 
much larger than those seen in numerical simulations, indicating that additional physical input 
to the fat-center-vortices model is required. At large distances, screening occurs for zero-triality 
representations; for the representations with non-zero triality the string tension equals that of the 
fundamental representation. Some rather "unphysical" flux distributions can lead to violations of 
f"^ ■ Casimir scaling at intermediate distances and violations of the expected ordering of representations 

f*^ ) at large distances. 
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I. INTRODUCTION 



> 
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A variety of models have been introduced in an attempt to explain quark confinement. The center vortices model 
introduced in the late 1970's by 't Hooft is one of those attempts. A center vortex is a topological line-like (in 
D=3 dimensions) or surface-like (in D=4 dimensions) field configuration. The vortex carries magnetic flux quantized 

■ in terms of elements of the center of the group. The fluxes form narrow tubes with constant energy per unit length 
| (surface). In order for the vortex to have finite energy per unit length, the gauge potential at large transverse 

distances must be a pure gauge. However, the gauge transformation which produces that potential is non-trivial. It 
is discontinuous by an element of the gauge center. It is the non-trivial nature of the gauge transformation which 
forces the vortex core to have non-zero energy and makes the vortex topologically stable. Faber, Greensite and 
Olejnik j|] introduced fat-center-vortices to obtain confinement of both fundamental and higher representation static 
, sources. According to the fat-center-vortices model, the vacuum is a condensate of vortices of some finite thickness. 

■ Confinement is produced by the independent fluctuations of the vortices piercing each unit area of a Wilson loop. 
Faber et al. explicitly worked out the model for SU(2) and gave results for a particular flux distribution within 

the vortices. Here, I work in SU(3) and investigate a wide variety of vortex flux distributions. I study the existence 
of Casimir scaling at intermediate distances and the pattern of screening at large distances. The model confirms 
Oh, Casimir scaling qualitatively at intermediate distances, however, in general the model does not agree with Casimir 
scaling quantitatively, in contrast to the precise agreement recently found by numerical simulations [@j8|. The lack 
of convexity of the potentials predicted by the model is also discussed. For completeness, I first briefly explain their 
model and then apply it to SU(3), using it to study the potentials between static quarks for the fundamental and 
several other representations. 



II. THE MODEL OF FABER, GREENSITE AND OLEJNIK 



In the fundamental representation of SU(N), a center vortex linked to a Wilson loop has the effect of multiplying 
the Wilson loop by the gauge group center, 

W(C) -> exp 2 ^ W{C) n = 1,2,..., N- 1. (2.1) 

Based on the vortex theory, the area law for a Wilson loop is due to the quantum fluctuations in the number 
of center vortices linking the loop. Adjoint Wilson loops are not affected by center vortices unless the vortex core 
overlaps the perimeter of the loops. If the vortex thickness is large enough, the fat-center-vortices model can explain 
confinement and the Casimir scaling of higher representation string tensions. The average Wilson loop predicted by 
this model has the following form: 

N-l 

< W(C) >= - E M 1 - ReGr{c%(x v )})}, (2.2) 

x v n—1 
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where x v is the location of the center of the vortex, A is the area of the loop C, and Q r is defined as: 

g r [a] = ^-Trexp[ia.H}. (2.3) 

d r is the dimension of representation r, and {Hi} is the subset of the generators needed to generate all elements of 
the center of the group. (For SU(3), As is sufficient.) The parameter / represents the probability that any given unit 
area is "pierced" by a vortex; i.e., a line running through the center of the vortex tube intersects the area. 

The parameter a c (x) is determined by the fraction of the vortex flux that is enclosed by the Wilson loop. Therefore 
a c (x) depends on the profile of the vortex as well as the shape of the loop and the position x v of the center of the 
vortex relative to the perimeter. For SU(3) a c (x) is equal to 47r/-\/3, if the flux is entirely inside the minimal area of 
the loop and it is zero if the flux is entirely outside the minimal area of the loop. 

III. GENERAL FEATURES OF CASIMIR SCALING AND SCREENING 



Numerical simulations |l|-|8|] have shown that the potentials in SU(3) at zero temperature are linear at intermediate 
distances and roughly proportional to the Casimir operator. The proportionality of potential to the Casimir operator 
is known as "Casimir scaling." Assuming such scaling, the potentials of each representation to that of fundamental 
representation would be 2.5, 2.25, 4.5, 7, 4 and 6 for representations 6, 8, 10, 15-symmetric, 15- antisymmetric and 
27, respectively. The Casimir scaling regime is expected to exist for intermediate distances, perhaps extending from 
the onset of confinement to the onset of screening. 

Screening can be understood as follows: Each representation can be labeled by the ordered pair (n, m), with n and 
m the original number of 3 and 3 which participated in constructing the representation. Triality is defined as (n-m) 
mod 3. Screening occurs for representations with zero triality: 8 = (1,1), 10 = (3,0), and 27 = (2,2). For these 
representations, as the distance between the two adjoint sources increases, the potential energy of the flux tube rises. 
A pair of gluons pops of vacuum when this energy is equal or greater than the twice of glue-lump mass. (A glue-lump 
is the ground state hadron with a gluon field around a static adjoint source.) At large distances, the static sources 
combine with the adjoint(8) charges (dynamic gluons) popped out of the vacuum and produce singlets which screen. 
Therefore the potential between static sources is no longer R dependent. Static sources in representations 10 and 27 
transform into the 8 first by combining with a dynamic gluon, and then the 8 transforms into the singlet by combining 
with a second gluon. 

8®8 = 27efoeioe8©i, (3.1) 

10® 8 = 8® 100 27 35, (3.2) 

27 <g> 8 = 64 © 27 © 27 © 35 © 3~5 © 10 © 10 © 8. (3.3) 

Therefore, we expect the potential in representations 10 and 27 to screen only at higher energy than in representation 
8. Static sources in representations with non-zero triality, 6 = (2, 0), 15 s = (4, 0) and 15 a = (2, 1), transform into the 
lowest order representation (3 and 3) by binding to the gluonic 8's which are popped out of the vacuum: 

6®8 = 3©6©15ffi24, (3.4) 
15„ © 8 = 42 © 24 © 15 Q © 15 a © 6 © 3 © 15 s , (3.5) 

15 s ®8 = 48ffi42ffil5 s ffil5 a . (3.6) 

15-symmetric changes to 15- antisymmetric first, so it needs to interact with the 8's (popped from the vacuum) twice 
to transform to 3. Screening does not occur for representations with non-zero triality, since there is no way to get 
a zero triality representation by crossing a non-zero one with any number of 8's. As a result, the slope of the linear 
potentials of the representations with non-zero triality changes to the slope of the fundamental one, and a universal 
string tension is observed for large R. We expect the representation 15-symmetric to require a larger value of R to 
approach the fundamental slope than representations 6 or 15- antisymmetric because two pairs of 8's must be popped 
from the vacuum in the 15-symmetric case. 
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IV. APPLYING THE FAT-CENTER- VORTICES MODEL TO SU{3) 



To find the potential V r (R) in SU(3), first I need to find Hi in Eqn. 2.3 for each representation: 3, 6, 8, 10, 15- 

2 



symmetric, 15- antisymmetric, and 27. For the fundamental representation, H\ — Tg = 4^; where As is the diagonal 



Gell-Mann matrix. 

I obtain Tg of other representations by using the tensor method. Define {X' r ;i = 1, ...,c? r }, which are the basis 
vectors for the space on which the representation act. The corresponding generators are obtained from ||: 

Pa 2 ]ix,jy = [T a 1 ]ijS xy + 5j j \T a 2 ]xy (4-1) 

T' a s are the group generators for representations D\, D2, D\ ® T>i- The elements of Tg can be found by: 

d r 

Y^c^x-. (4.2) 

To study potentials for SU(3), one needs to define an appropriate form of the function a c (x) in Eqn. 2.2. To 



understand more about the effect of the vortex profile on potentials and Casimir scaling, I assume a density of flux 
p(r) in an axially symmetric vortex core, where r is the radial distance from the vortex center. Let p(r) — for r > a 
so the vortex has a sharp boundary at r — 0. Now let (3(x v ) denote the amount of the flux of the vortex contained in 
the region x > 0. Thus, (3(x v ) = for x v < —a, and (3{x v ) — ^= for x v > a. For —a < x v < a, f3{x v ) is determined by 

the integral of p(r) over the fraction of the vortex in the region x > 0. Finally, let the Wilson loop have sides x = 
and x = R. Then a c {x v ), the fraction of flux within the loop, is given by: 

ocr{x v ) = (3{x v ) - /3(x v - R), (4.3) 

A simple choice for p(r) that I have tried is a uniform distribution p(r) = po for r < a. Another possibility with a 
smoother edge at r = a, is: 

p(r) = po exp[ b - }, (4.4) 

(a !) 

where b is an adjustable constant and po is fixed by the requirement that the total flux is Fig. [l] shows potentials 
obtained from this flux distribution with b = 0.1, a = 20 and / = 0.1. From the plot, it can be seen that, for each 
representation, there exists a region in which the potential is approximately linear and qualitatively in agreement with 
Casimir scaling. Screening occurs for representations 8, 10 and 27 while the slope of the potentials for representations 
6, 15-symmetric and 15-antisymmetric changes to the slope of the fundamental representation. Note the non-convexity 
near R = 8 to R — 20 for all representations, and especially representations 15-symmetric and 27. Even though the 
fat-center-vortices model predicts some of the expected behavior of the potential between static quarks, it has some 
limitations. In particular, it violates the fact that the potential should be always a convex function of distance |^o| . 

Qualitative agreement with Casimir scaling is observed for all the axially symmetric distributions I tried. This is 
true even if one defines the density to be zero everywhere except on the outer boundary (r — a) of the vortex. Fig. 
H plots potentials for this distribution with a = 20 and / = 0.1. Fig. || shows potentials for the maximally non- 
axially symmetric core where the flux is zero everywhere except at the two points where the vortex first enters and 
exits the Wilson loop. A linear regime still exists at intermediate distances but qualitative agreement with Casimir 
scaling is lost: slope of potentials with larger Casimir operators have smaller string tensions. For example, the string 
tension for representations 6 and 8 are larger than the ones for representations 15-antisymmetric and 27. In this 
case the order of potentials at long distances changes as well. For example, the potential for representation 27 is less 
than the potential for representation 8 in the screened regime. It still remains true, however, that the zero triality 
representations screen, and non-zero triality ones approach the fundamental string tension at long distances. However 
the non-convexity of potentials is almost gone. Note that this flux distribution is probably unphysical. We expect 
the flux in the lowest energy vortex to be axially symmetric. The fact that potentials obtained from physical flux 
distributions agree qualitatively with Casimir scaling, is the strength of the vortex model. However recent numerical 
simulation results B show that potentials are quantitatively in agreement with Casimir scaling with an accuracy of 5 
percent, the feature that is lost in fat-center-vortices results. For example, from Fig. refPdens, ratios of potentials for 
representations 6, 8, 10, 15-symmetric, 15-antisymmetric, and 27 to the fundamental representation are 2, 1.82, 3, 5.4, 
2.7, and 3.4, respectively whereas the same ratios for Casimir scaling would be 2.5, 2.25, 4.5, 7, 4, and 6, respectively. 

One can get similar results for potentials at intermediate and long distances using one dimensional functions for 
a(x) similar to those of ref. 0. I have tried several functions. An example is: 
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(3(x v ) 




v > a 

V ^ 

-a < x v < a. 



- a 



(4.5) 



&r(x v ) is then given by Eqn. 4.3 
is needed to find 



Note that in this case no density distribution is defined and therefore no integration 
For moderate b this results in a flux profile similar to the one used in ref. || and leads to 
qualitative Casimir scaling. However for b large enough this gives similar plot to Fig. ^ and violates Casimir scaling. 
Thus integrating a circularly symmetric distribution inside the Wilson loop is more physical than just picking a 
function (3(x) at random. An arbitrary one dimensional distribution is not necessary consistent with axial symmetry, 
and in the above example, large b certainly does not correspond to any axially symmetric distribution. 



V. CONCLUSION 



By applying the fat-center- vortices model to SU (3) and using presumably physical axially symmetric density dis- 
tributions for the vortex, I showed that for several representations there exists a region at intermediate distances in 
which the static potential is linear and qualitatively in agreement with Casimir scaling. This is also in agreement 
with the observation in SU(3) simulations of a linear potential in proportion to Casimir ratio of the representation 
0-0. However, the Casimir proportionality is dependent on the flux distribution in the vortex and it is possible to 
lose this feature by changing the distribution function to a non-axially symmetric distribution. At large distances, 
zero-triality representations will be screened and the potentials for non-zero triality representations parallel the one 
for the fundamental representation. Some of the expected features of the screening pattern are also lost for non-axially 
symmetric distributions. The conclusion is that Casimir scaling and the pattern of screening depend on the detailed 
vortex structure and are not simple kinematic consequence of the fat center vortex picture. However it is also clear 
that these properties are rather robust and are likely to survive with physically realistic vortices. On the other, po- 
tentials are not quantitatively in agreement with Casimir scaling as predicted by recent numerical simulations. This 
suggests that the fat-center- vortices model needs further refinement if it is to remain viable. In particular, one may 
need an appropriate physical flux distribution of vortex sizes. Further numerical studies of these issues is in progress. 
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FIG. 1. Potentials for different representations using density distribution proportional to p exp( ^ |T , | 2 ) for —a < x < a. 

The fundamental representation is shown by the letter "f". At intermediate distances (3 < R < 5), the ratios of the string 
tension of representations 8, 6, 15-antisymmetric, 10, 27 and 15-symmetric to the one of the fundamental representation are 
2.02, 2.21, 3.1, 3.4, 3.8 and 5.6, respectively. This is qualitatively in agreement with the Casimir ratios which are 2.25, 2.5, 4, 
4.5, 6 and 7, respectively. 
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FIG. 3. Static sources potential using a density distribution proportional to delta function which is zero every where except 
at the two points where the vortex first enters and exits the Wilson loop. 
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